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SUMMARY 

The purpose of tht. paper la to present stress-intensity factors for a 
wide range of nearly seml-elllptlcal surface cracks In pipes and rods. The 
surface cracks were oriented on a plane normal to the axis of pipes or rods. 

The configurations were subjected to either remote tension or bending loads. 

For pipes, the ratio of crack depth to crack length (a/c) ranged from 0.6 
to 1; the ratio of crack depth to vail thickness (a/t) ranged from 0.2 to 0.8; 
and the ratio of Internal radius to wall thickness (R/t) ranged from 1 to 10. 
For rods, the ratio of crack depth to crack length also ranged from 0.6 to 1; 
and the ratio of crack depth to rod diameter (a/D) ranged from 0.05 to 0.35. 
These particular crack configurations were chosen to cover the range of crack 
shapes (a/c) that have been observed In experiments conducted on pipes and 
rods under tension and bending fatigue loads. The stress-intensity factors 
were calculated by a three-dimensional finite-element method. The finite 
element models employed singularity elements along the crack front and linear- 
strain elements elsewhere. The models had about 6500 degrees of freedom. The 
stress-intensity factors were evaluated using a nodal-force method. 

The present results were co, spared with other analytical and experimental 
results tor some of the crack configurations. The re.ults generally agreed 
within 10 percent. 


♦Senior Scientist, Analytical Services & Materials, Inc, Tabb, VA 
**Senior Scientist, Materials Division. 
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Th e present results will be useful to predict crack-growth rates and 
fracture strengths. to design structural components, and to establish Inspec- 

tion intervals. 

INTRODUCTION 

Surface cracks can occur In many structural components. Circumferential 
surface cracks can cause premature failure of landing gear of aircraft 
piping, bolts, pins, and reinforcements which employ cylindrical shaped com- 
ponents. .curate stress analyses of these surface-cracked components are 
needed for reliable prediction of crack-growth rates and fractu.e strong • 
However, because of the complexities of such problems, exact solutions are 

available. - 

some investigators have used experimental or approximate analytica 

methods to obtain stress-intensity factors for surface cracks In rods under 

tension and bending loads. The stress-intensity factors for a circumferential 

„ al mi us in K the James-Anderson procedure 12]. 
experimentally by Wilhem, et al. U1 usin & 

. intersect the outer surface of the 
Their results show that the surface cracks 

. ot a i. [3], (boundary-integral 

rod at nearly right angles. Athanasslaols , 

, 141 (finite-element method) used analytical methods 

method) and Netu, et al. (41, 

a p p nrfi for circumferential surface cracks ot 
to obtain stress-intensity factors for circum 

various shapes in rods. Their surface-crack shapes were chosen to agree « 
Chair experimental observations. In some o, the experiments of Reference 4, 
however, the surface cracks did not Intersect the outer surface at right 
angles. Residual stresses at the outer surface could have altered the crac 
Shape near the free surface. Trantlna. et al. .5. conducted an elastic and an 
elastlc-plast ic finite-element analysis of a small surface crack in a rod 
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under tension loading. They assumed that cracks intersected the outer surface 
at right angles. 

Much less research has been conducted on circumferential surface cracks 
in hollow cylinders or pipes than on rods. Delale and Erdogan [6] and German, 
et al. [7] obtained stress-intensity factors for interior and exterior circum- 
ferential surface cracks using the line-spring model. Recently, Forman and 
Shivakumar [8] studied the fatigue-crack growth behavior of circumferential 
surface cracks in rods and pipes. Most of their experimental results did in- 
dicate that surface cracks intersect the free surfaces at about right angles. 

This paper presents stress-intensity factors calculated with a three- 
dimensional finite-element analysis for a wide range of nearly semi -elliptical 
surface cracks in pipes and rods. The surface cracks were oriented on a plane 
normal to the axis of pipes or rods, as shown in Figure 1. The crack 
configurations were assumed to be such that the crack fronts intersect the 
free surface at right angles. The pipe and rod were subjected to either 
remote tension or bending loads. For pipes, the ratio of crack depth to crack 
length ranged from 0.6 to 1; the ratio of crack deptn to wall thickness ranged 
from 0.2 to 0.8; and the ratio of internal radius to wall thickness was 1 to 
10. For rods, the ratio of crack depth to crack length also ranged from 0.6 
to 1 ; and the crack configurations were chosen to cover the range of crack 
shapes that have been observed in experiments conducted on pipes and rods 
under tension and bending fatigue loads. The stress-intensity factors were 
calculated by a nodal-force method (9-11]. The present results were compared 
with other analytical and experimental results from the literature for some of 
the crack configurations. 
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SYMBOLS 


depth of surface crack 

half-arc length of surface crack 

i outer diameter of pipe and rod 

' stress-intensity houndary-correction factor 

boundary-correction factor at maximum depth point on surface 
crack (<|> - ir/2) 

boundary-correction factor at intersection point of crack and 
outer surface 

1 half-length of pipe and rod 

l stress-intensity factor (mode I) 

C A stress-intensity factor at maximum depth point on surface 

crack (<J» = tr/2) 

} shape factor for elliptical crack 

R. internal radius of pipe 

remote bending stress on outer fiber 
S t remote uniform-tension stress 

t wall thickneas of pipe 

x,y Cartesian coordinate system 

v Poisson's ratio 

<(> parametric angle of ellipse 

Superscript 

• primes denote quantities associated with surface crack in flat plate 

THREE-DIMENSIONAL FINITE-ELEMENT ANALYSIS 
A three-dimensional finite-element analysis v*s used to calculate the 
mode I stress-intensity factor variations along the crack front for a circum- 
ferential surface crack in the pipe and rod shown in Figures 1(a) and 1(b), 
respectively. The pipe and rod was subjected to either remote tension or 



bending loads. In this analysis, Poisson’s ratio (v) was assumed to be 0.3. 

The shapes of the surface cracks were nearly but not exactly semi-elliptical. 
These crack shapes were generated using a conformal transformation as 
described in tne appendix. 

Figures 2(a) and 2(b) show typical finite-element models for a circumfer- 
ential surface crack in a pipe and rod, respectively. The finite-element 
models employ singularity elements along the crack front and linear-strain 
elements elsewhere. The models had about 6500 degrees of freedom. Stress- 
intensity factors were evaluated from a nodal-force method. Details of the 
formulation of these types of elements and of the nodal-force method are given 
in References 9-11 and are not repeated here. Details on the development of 
the finite-element models are given in the appendix. 

Loading 

Two types of loads were applied to the finite-element models of the 
surface-cracked pipe and rod: remote uniform-tension and remote bending. The 

remote uniform-tension stress is S t and the remote outer-fiber bending 
stress is S b . The bending stress S b , in Figure 3, is calculated at the 
origin of the surface crack (x = y = 0 in Fig. 4) without the crack being 

present. 

Stress-Intensity Factor 

The tension and bending loads only cause mode I deformations. The mode I 
stress-intensity factor K for any point along the surface-crack was taken to 

be 


K = 



( 1 ) 
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where the subscript l denotes either tension load (i - t) or bending load 
(i * b), and Q, the shape factor for an ellipse, is given by the square of 
the complete elliptic integral of the second kind. The half-length of the 
pipe and rod, h, was chosen large enough to have a negligible effect on 
stress intensity (h/D > 10). Values for F, the boundary-correction factor, 
were claculated along the crack front for various combinations of parameters 
(a/t, a/c, R/t, and $ for a crack in a pipe; and a/c, a/D, and <|> for a 
crack in a rod). The crack dimensions and parametric angle, <J>, are defined 
in Figure 4. The range of crack shapes (a/c) and of sizes (a/t or a/D) 
analyzed are shown in Figures 5 and 6 for the pipe and rod, respectively. 

The empirical expressions for Q used in this paper were developed by Rawe 
(see Ref. 11) and are 

Q = 1 + 1.464(a/c) 1 *^ 5 for a/c < 1 ( 2a ) 

Q = 1 + 1.464(c/a) 1 *^ 5 for a/c > 1 i 2b ) 

RESULTS AND DISCUSSION 

In the following sections, stress-intensity factors for various shape 
surface cracks in pipes and rods subjected to tension and bending loads are 
presented. Tables 1-3 give the normalized stress-intensity factors, K/S/ira/Q, 
at the maximum depth point (A) and at the point which the crack intersects the 
free surface (B). Figures 7-12 show the variation in normalized stress- 
intensity factors, as a function of the parametric angle (2$/ir) for various 
crack shapes (a/c), crack size (a/t or a/D), and radius of the pipe (R/t). 


Pipes Under Tension Loads 

Figure 7 shows the normalized stress-intensity factors as a function of 
the parametric angle (+) for a pipe (R/t - 2) subjected to remote tension with 
a semi-circular crack (a/c = 1) for various values of a/t. For this crack 
shape, the maximum normalized stress-intensity factor occurred at the point 
where the crack meets the free surface (point B). For all of the a/c ratios 
considered (from 1 to 0.6), larger values of a/t always gave larger normalized 
stress-intensity factors. 

For an a/c ratio of 0.8 and a/t less than or equal to 0.5, the 
normalized stress-intensity factors at the deepest point and at the free 
surface are nearly the same (see Table 1). On the other hand, for an a/c 
ratio of 0.6, the maximum normalized stress-intensity factors occurred at the 

point of maximum depth (<(> = tt/2). 

Figure 8 shows the normalized stress-intensity factors for a surface 
crack with a/c = 0.8 and a/t =0.5 in pipes with R/t ratios of 1 and 
10. (The results for R/t = 2 and 4 lie in between those for R/t = 1 and 
10 and are not shown for clarity.) For this configuration, the effect of 
varying R/t is insignificant. However, for cracks with a/c =0.6 and 
a/t = 0.8, R/t has a significant effect on the normalized stress-intensity 
factors, as shown in Figure 9. The values at the deepest point (* = */2) are 
affected more than those at the free surface. Figure 9 shows that lower 
R/t values gave higher stress-intensity factors. However, the differences 
between the stress-intensity factors are less for larger values of R/t. Thus 
the effect of R/t diminishes for pipes with larger R/t values. 

In summary, the effect of the curvature of the pipe (R/t) is to elevate 
the stress-intensity factors compared to a flat plate (R/t = »). The effect 
is more pronounced at the deepest point than at the free surface point. 


Pipes Under Bending Loads 

The normalized stresa-intensity factors as a function of the parametric 
angle ($) are shown in Figure 10 for a pipe (R/t - 2) subjected to remote 
bending with a semi-circular crack for various values of a/t. For this crack 
shape, the deeper cracks (larger a/t ratios) produced larger normalized 
stress-intensity factors where the crack meets the free surface but smaller 
values at the maximum depth point (<|> * ir/2). 

Rods Under Tension Loads 

Figure 11 shows the normalized stress-intensity factors as a function of 
the parametric angle (*) for a rod with various shape surface cracks with 
a/D - 0.2. When a/c was unity, the maximum normalized stress-intensity 
factor occurred at the free surface. When a/c was equal to 0.6, the maximum 
was at the deepest point. For surface cracks with an a/c ratio of 0.8, 
however, the normalized stress-intensity factors are nearly constant, much 

like the pipe. 

Rods Under Bending Loads 

The normalized stress-intensity factors as a function of the parametric 
angle (<(.) for a rod with various shape surface cracks with a/D - 0.2 are 
shown in Figure 12. The maximum normalized stress-intensity factor occurred 
at the free surface for a/c - 1 and a/c » 0.8. For surface cracks with 
a/c * 0.6, however, the normalized stress-intensity factors all along the 
crack front are nearly constant. 

Comparisons With Other Solutions 

The comparison of the present results with those from the literature are 
difficult because, at least, three definitions of crack shapes have been used 
All definitions differ in how the crack length c is measured. In this 
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paper, c is measured as the arc length, as shown in Figure 4. In some 
reports, c is measured as the horizontal projection of point B on the 
x-axis, while in other reports, c is defined as the intersection point of an 
ellipse with the x-axis* In the latter case, the crack front will not inter- 
sect the free suface at a right angle. The crack front shape, in the latter 
case, therefore, will be different from that used in this paper. In view of 
these difficulties, only a few comparisons can be made. 

As previously mentioned, stress— intensity factor analyses of circumfer- 
ential surface cracks in pipes have received very little attention in the 
literature. Delale and Erdogan [6] obtained stress-intensity factors for 
interior and exterior circumferential surface cracks and German, et al. [7] 
obtained stress-intensity factors for interior circumferential surface cracks 
by using the line-spring model. Most of the external circumferential surface 
crack configurations presented in Reference 6 were vastly different from the 
configurations presented in this paper. However, one configuration with 
a/c = 0.775, a/t = 0.8 and R/t = 5.374 falls within the range of param- 
eters considered in this paper. (In the notation of Ref. 6, this configura- 
tion has L Q /h = 0.8, a/h = 1 and A 2 = 0.75.) The normalized stress- 
instensity factor at the deepest point (<p — it /2) of the crack was computed 
from the results of Reference 6 as 1.276. Interpolating the present results 
in Table 1, the normalized stress-intensity factor F A for this configuration 
was found to be 1.161. The result from the line-spring model, reference 6, is 
aoout 10 percent higher than the present result. 

The rod configurations with surface cracks subjected to remote tension 
have received more attention in the literature than the pipe configurations. 
These configurations were analyzed by Wilhem, et al. [1], Athanassiadis, 
et al. [ 3 J , and Nezu, et al. [4]. The present results are compared with the 
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results from References 1 and 3. Comparisons with the results from Nezu, 
et al. [4] could not be made because the crack shapes analyzed in Reference 4 
and those In the present analysis (see Pig. 6) were very different. 

Figure 13 compares the normalized stress-intensity factors at the free 
surface <P B ) and the maximum depth point (P A ) for a surface crack with 
a/c = 0.6 from the present finite-element analysis with those from a Boundary 
Integral Equation (B1E) method [3]. The results of Athanassiadis, et al. 13] 
were interpolated and plotted in Figure 13 as open symbols. The present 
results are shown by solid symbols. The normalized stress-intensity factors 

obtained by the BIE method were 0 to 10 percent lower than the present 
results. 

Figure 14 compares the normalized stress-intensity factors at the maximum 
depth point for surface cracks with various shapes (a/c) and sizes (a/D) from 
the present analyses and from experimental results. Wilhem, et al. [1] 
obtained an experimental stress-intensity factor solution using the James- 
Anderson procedure [2]. These results are shown by the dashed curve in 
Figure 14. For the surface cracks in their tests, the a/c ratios varied from 
0.95 to 0.85. The present results (symbols) for a/c * 0.8 and 1.0 bound the 
experimental results for a/D < 0.25 and are a little below for a/D > 0.25. 

Bush [12] considered cracks with straight fronts (see insert in Fig. 14) 
in rods subjected to remote tension. He obtained stress-intensity factors 
from experimental compliance for these straight through cracks of various 
depths. His results are shown in Figure 14 by a solid curve. For a surface 
crack with an a/c ratio of 0.6 and an a/D ratio of 0.35 (see Fig. 6(c)), 
the crack configuration is very nearly the same as that for a crack with a 
straight front. For this configuration, the present results for a/c - 0.6 
are a little below (about 2 percent) the straight through crack results. The 


10 


present results need not necessarily agree with the experimental results 
because the crack shapes are not Identical, as noted previously. 

CONCLUDING REMARKS 

Stress-intensity factors for circumferential surface cracks In pipes and 
rods have been obtained by a three-dimensional finite-element analysis. The 
pipes and rods were subjected to either temote tension or remote bending 
loading. The surface cracks were nearly semi-elliptical and were oriented on 
a plane normal to the axis of pipes or rods. A wide range of crack shapes, 
crack sizes, and internal radius-to-wall thickness ratios have been con- 
sidered. For each of these crack configurations and loadings, the stress- 
intensity factors calculated by the finite element analysis are presented. 

Stress-intensity factors for surface cracks in a pipe were found to be 
insensitive to internal radius-to-wall thickness (R/t) ratios ranging from 
1 to 10, for crack depth-to-length (a/c) ratios ranging from 0.8 to 1.0 with 
crack depth-to-wall thickness (a/t) ratios less than 0.8. For a/c =0.6 
and a/t = 0.8, however, the stress-intensity factors showed significant 
variation with R/t. The effect of the curvature of the pipe (R/t) is to 
elevate the stress-intensity factors compared to a flat plate (R/t = «). This 
effect is more pronounced at the deepest point than at the free surface point. 

Stress-intensity factors for a surface crack in a rod were 0 to 10 per- 
cent higher than those calculated from a boundary-integral analysis. The 
stress-intensity factors agreed well with experimental results for surface 
cracks in rods and approached the experimental results for cracks with 
straight fronts. 

The stress-intensity factors obtained here should be useful in predicting 
fatigue-crack growth and fracture of surface cracks in cylinders and rods. 
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APPENDIX 


The purpose of this appendix Is to present the procedure used to develop 
finite-element models for surface cracks in pipe and rod configurations 
through a conformal transformation. 

A cylinder with a surface crack is shown in Figure 15(a). The stress- 
intensity factors for this configuration were evaluated from a nodal-force 
method [9]. In this method, the nodal forces normal to the crack plane (x,y 
plane) and ahead of the crack front are used. The nodal- force method also 
requires that these forces be evaluated at nodes which are very near the crack 
front and which lie on lines in the x,y plane that are normal to the crack 
front. Therefore, the finite-element model should be such that the normality 
at the crack front is maintained. This is achieved through a conformal trans- 
formation as follows. 

First, a finite-element model for a semi-elliptical surface crack with 
semi -minor and semi -major axes, a' and c'» respectively, in a plate of 
width w' and a thickness of t* (see Fig. 15(b)) is developed such that 


a’ 



R + t 
+ t - 


t) 


(3) 


R + t 


(A) 


and 


t' 


In 


R + t 
R 


(5) 


To obtain the desired configuration in Figure 15(a), a conformal trans 
formation 


x 


(R + t) e" y cos(w'/A - x') 
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( 6 ) 


< *> « • 




y ■ (R + t)[l - e~y' 8in(w*/A - x' ) 


(7) 


z 


z’ 


and 

w' * 2ir (8) 

is used. This transformation transforms every point in the x'.y'.z' system 
to a unique point in the x,y,z system and maintains normality. Because the 
finite-element model of a surface crack in a flat plate (Fig. 15(b)) has nodes 
along hyperbolas near the crack front [11] (and, hence, normality to the semi- 
elliptical crack front in the x',y' plane is assured), the conformal trans- 
formation gives nodes along curves in the x,y plane which are also normal to 
the crarW front in the pipe configuration (Fig. 15(a)). 

The finite-element models for the surface crack in the rod configuration 
were obtained from the pipe models by idealizing the inside core with finite 
elements (see Fig. 2). 
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Table 1.- Normalized stress-intensity factor, K/S /ira/Q, for surface crack 

in a pipe subjected to tension 'loads. 

A * maximum depth point 
B » free surface point 





































Table 2.- Normalized otress-intensity factor, K/S^/na/Q, surface crack 
In a pipe subjected to bending loads. 

A = maximum depth point 
B = free surface point 





























Table 3.- Normalized stress-intensity factor, K/S/ira/Q, for surface crack 
in a rod subjected to tension or bending loads. 

A =» maximum depth point 
B = free surface point 


Tension loads 


a/c = 0.8 



a 

a/c = 1.0 

a/c * 0.8 

a/c * 0.6 

D 

A 

B 

A 

B 

A 

B 

.05 

0.938 

1.129 

0.984 

1.029 

1.035 

0.907 

.125 

0.836 

1.114 

0.901 

1.019 

0.987 

0.903 

1.20 

0.749 

1.112 

0.830 

1.028 

0.985 

0.909 

1.275 

0.683 

1.109 

0.795 

1.040 

1.041 

0.924 

1.35 

0.629 

1.106 

0.782 

1.039 

1.056 

0.876 
























. coordinate system and dimensions 


for surface crack 
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Surface-crack shapes and sizes analyzed for 






Normalized stress-intensity factors along the front of a surface 
crack (a/c = 1) in a pipe under tension. 
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Fig. 15.- Conformal transformation of a surface crack in a plate to a surface crack 
in a pipe. 
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